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We consider a problem in quantum theory that can be formulated as an optimisation problem 
and present a global optimisation algorithm for solving it, the foundation of which relies in turn on 
a theorem from quantum theory. To wit, we consider the maximal output purity v q of a quantum 
channel as measured by Schatten (/-norms, for integer q. This quantity is of fundamental importance 
in the study of quantum channel capacities in quantum information theory. To calculate v q one has 
to solve a non-convex optimisation problem that typically exhibits local optima. We show that 
this particular problem can be approximated to arbitrary precision by an eigenvalue problem over 
a larger matrix space, thereby circumventing the problem of local optima. The mathematical proof 
behind this algorithm relies on the Quantum de Finetti theorem, which is a theorem used in the 
study of the foundations of quantum theory. 

We expect that the approach presented here can be generalised and will turn out to be applicable 
to a larger class of global optimisation problems. We also present some preliminary numerical 
results, showing that, at least for small problem sizes, the present approach is practically realisable. 



I. INTRODUCTION 

In Quantum Information Theory (QIT), noisy com- 
munication channels are modelled as completely positive 
trace-preserving maps between operator algebras. One 
of the most fundamental questions in QIT is the deter- 
mination of the capacity of these quantum channels to 
transmit classical information, the so-called classical ca- 
pacity of a quantum channel [8]. What makes the de- 
termination of this quantity so much more difficult than 
its purely classical counterpart is the existence of entan- 
glement in quantum physics. In fact, it has been shown 
that to obtain an optimal quantum channel decoder one 
has to perform entangled measurements over the channel 
output states. What is not known, however, is whether 
entanglement is also necessary to obtain an optimal en- 
coder. This is widely believed not to be the case, i.e. no 
benefit is expected in having entanglement between the 
single-letter states. To prove this, it is necessary to show 
that the single-letter classical capacity (a.k.a. the Holcvo 
capacity x) of a quantum channel is additive. 

The Holevo capacity of a channel $ is defined by 

x m = supscyv^o*)) -yv ( s($( ft )), 

7T,p 

where the sup runs over all probability distributions {iti} 
and collections of states {pi}, and S is the von Neumann 
entropy S(p) = — Tr[/>logp]. The additivity of this ca- 
pacity would be the property x($i®3>2) = x(®i)+x(&2), 
i.e. there would be no benefit in sending entangled states 
through the tensor product channel. 

The expression for the Holevo capacity looks quite 
complicated, and it has resisted any attempt so far at 
proving its additivity. For that reason, a simpler channel 
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property, called the minimal output entropy (MOE) [1], 
has been introduced, in the hope that proving additivity 
for the MOE would be simpler and yet shed some light 
on the additivity problem for the Holevo capacity. The 
minimal output entropy v$ of a channel <I> is defined as 

v s {§) = min5'(4>(p)), 

p 

where the minimum runs over all density matrices (i.e. 
unit trace positive semidefmite matrices). The entropy is 
minimal for those output states that are closest to being 
pure (rank 1). It also makes sense to consider the related 
quantities 

v q {§) = maxp(p)||g, 
p 

where purity is measured by the Schatten g-norm. These 
quantities are called the maximal output purities (MOP) 
and they are directly related to the MOE: 

1 - !/*($) 

!/<?($) =lim q —— . 

V ' ?U q-l 

Furthermore, if the MOP is multiplicative, Vq(Qi<S>&2) — 
v q {<&i)v q {<&2), or if it is at least for values of q close to 
1, then the additivity of the MOE follows as a direct 
consequence. 

Quite surprisingly, it has been shown recently that the 
additivity of the MOE is actually equivalent with the 
additivity of the Holevo capacity [16], in spite of the ap- 
parent greater simplicity of the expression for the MOE. 
What is more, additivity of the MOE is also equivalent 
with the additivity of the Entanglement of Formation 
[3, 16], another open problem in QIT whose resolution 
is eagerly awaited. Additivity of the MOE and multi- 
plicativity of the MOP have been proven in specific case 
[2, 11, 12, 13]. Unfortunately, there is a counterexample 
to the multiplicativity of the MOP for values of q > 4.79 
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[19], quenching the hope for a relatively simple proof of 
additivity of the MOE. Nevertheless, MOP might still be 
multiplicative for smaller values of q. Intuition has it that 
it might hold for q between 1 and 2. Since the Schatten 
(/-norms are easier to work with for integer values of q, 
it is reasonable to first try and prove multiplicativity of 
i>2 [14]. If V2 would indeed be multiplicative this would 
increase our belief in the above intuition. 

In the following, we study the quantity v q for integer 
values of q, and the value 2 in particular. The actual cal- 
culation of v q is an optimisation problem that is not very 
well-behaved, in the sense that it exhibits local optima. 
Indeed, it involves the maximisation of a convex func- 
tion over a convex set, and the only reasonable statement 
that convex analysis can make is that the maxima will be 
obtained in extremal points, i.e. pure states. Our main 
technical result is that the optimisation problem defining 
v q can for integer q be approximated by an eigenvalue 
problem over a larger Hilbert space. This approximation 
yields an upper bound and the error of the approximation 
goes to zero when the dimension of the enlarged Hilbert 
space goes to infinity. 

II. MAIN RESULTS 

Consider a c?-dimensional Hilbert space H, with its as- 
sociated state space S(H). We will denote the set of 
bounded Hermitian operators over Ti by B H (TL). Thus 
S(Ti) is the subset of B H (TL) containing all unit trace 
positive semidefinite operators. 

In the following, we will consider the quantity Vq(^) = 
max pe 5(-ft) Tr[(<f>(/?)) 9 ], for integer q. The maximand can 
be written as Tr[$(p)$(p) . . . <&(p)], with q factors. To 
every completely positive map $ we can associate a posi- 
tive semidefinite block matrix <fr, called the Choi matrix, 
such that &(p)ij = J2k i^kiPkl- Hence, the maximand 
can be brought in the general form Tr[Ap® 9 ], where A 
is a matrix in M^M^.-.M^C)...)) (with q occurrences 
of Md) depending on $. This matrix A has compos- 
ite row and column indices which we will denote by 
(i) := («i,«2, ■ ■ ■ ,i q ) and (J). We will think of A as be- 
ing an operator over the tensor product 7i® 9 , which is a 
Hilbert space consisting of q copies of H. Explicitly, we 
have 

A (i),U) =^[**iJi---*^J- (!) 

For q = 2, A is Hermitian, for higher q this is generally 
not true. 



A. Symmetry 

Consider permutations of n copies of Ti, 
where S n is the symmetric group of order n. If, for ev- 
ery permutation ir e S„, a matrix A over 'H® ri obeys 

A (il,...,i g ),{jl,-J<l) = AV(1).->V(,)).0V(1).-.JV(9))' tliei1 ^ 



matrix A is symmetric (not to be confused with transpo- 
sition symmetry) . To denote the action of a permutation 
7r G S n on a composite index, we will use the abbrevi- 
ation ir(i) := (vm, • • • ) i-K(n))- Let P-rr be the permu- 
tation matrix that permutes the indices according to ir, 
i.e. {P v x)u-s = XnU)- Thus A is symmetric if and only if 
W G S n , P\APk = A. The set of all symmetric Hermi- 
tian operators over H® n will be denoted by B HS {H® n ). 
This set is obviously a subspace of B H (H® n ). The set 
of all symmetric states over will be denoted by 

S s (H® n ). The state p® q is a simple example of a sym- 
metric state. The linear map P„ that projects all opera- 
tors in B{H® n ) to B s {H® n ) is given by 

We call P n {A) the symmetric part of A. 

If, as in our case, the matrix A is defined by the re- 
lation (1) where $ is a completely positive map, then A 
is in general only Hermitian and symmetric for q = 2. 
However, P q (A) is Hermitian for all values of q (and, of 
course, symmetric), as follows easily from (1). 

We will also need to consider the symmetry properties 
of vectors. Vectors tjj G 'H® n are totally symmetric if, 
for all 7r G S n , P-nip = V 1 ! or V'(i) — VV(i)- The totally 
symmetric vectors form a subspace, which we will denote 
by S(H <8n ), and which has dimension S(d, n) := CJ+f" 1 . 
We denote by S„ the projector on S(H® n ), that is 

s = — V p 

Finally, P n denotes the matrix whose columns span 
S(H® n ) and arc normalised such that P r \P n = 1. 

B. Main Theorem 

Let 1 be the identity matrix over Ti. Denote the maxi- 
mal eigenvalue of a matrix X and its corresponding eigen- 
vector by X max (X) and ip max (X). 

Theorem 1 For any q, n G N, and for any A G B{H m ) 
with Hermitian symmetric part (P q {A) G B H (H® q )), 
there exists a non-increasing sequence (fi n (A)) n , with 

p n (A) :=A max (P ?+ „(A®l®")), 

converging to 

lim Hn{A) = maxTrUp® 9 ]. 

n^oo p 

Moreover, the optimal p is given by 

p pt = lim Tr q+n -i \4> n ) (4>n\ 

n— >oo 

= ^(P,+„(4®i® n )). 
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Note that the matrix A need not be symmetric; one of 
the things the Theorem tells us is that only its symmetric 
part is of relevance. 

The main ingredient of the proof of this Theorem is 
the so-called Quantum de Finetti Theorem (QdF) [9, 17], 
which is a non-commutative analog of de Finetti's The- 
orem from probability theory. The QdF theorem char- 
acterises the so-called exchangeable sequences of states 
(p^)ni which are sequences of symmetric states over 
H® n such that Vn,77i £ N : p^ = Tr TO [p( n+m )]. Here 
Tr m denotes the partial trace over m copies of TL. The 
QdF theorem states that a sequence is exchangeable if 
and only if there exists a positive measure dp(p) over 
S(H) such that for all n £ N: 



» 



/ P® n dp( P ). 

JS(H) 



The QdF theorem has been used in the study of the foun- 
dations of quantum mechanics [4, 10], in mathematical 
physics [7], and recently also in QIT [6]. 



= lim max{Tr[A Tr„[p]] : p £ S s \H® {q+n) )} 

n^>oo p 

= lim max{Tr[(A ® 1®") p] : p £ S s (H®<> q+n) )}. 

ro^oo p 

This maximisation over symmetric states can be replaced 
by a maximisation over all states, provided they are pro- 
jected first onto the symmetric subspace: 

maxTr^p® 9 ] 
p 

= lim max{Trp®l®")P„ +n (p)] : p £ S(H® {q+n) )}. 

n — >oo p 

The cyclicity of the trace of a matrix product implies that 
the projection can equally well be applied to the factor 
A <g> 1®", yielding 

maxTr[Ap®«] 
p 

= lim max{Tr[P„ + „(^ ® l® n ) p] : p £ S(H® {q+n) )} 

n— *oo p 

= lim \ max (P q+n (A®t® n )). (2) 



C. Proof of Theorem 1. 

We are concerned with the maximisation of Tr[^4p® 9 ] 
over all density matrices p. We first turn this into the 
more general optimisation problem 



max {Tr[^4 

dp(p)>o 



p® q d»(p)] : J dp(p) = 1}, 



which obviously achieves the same maximum value; the 
optimal measure dp will be the Dirac measure peaking 
at the optimal p. We now claim that we can replace 
the latter maximisation by a maximisation of Trf-Ap^] 
over all p(«) that are g-th element in some exchangeable 
sequence of states. Indeed, any state of the form p® q 
obviously forms part of an exchangeable sequence. Con- 
versely, by the QdF theorem, the g-th element in any 
exchangeable sequence must be of the form J p® q dp(p) 
for some positive measure dp,{p). 

From the definition of exchangeable sequence, we infer 
that p( ? ) must be the partial trace of a symmetric state 
p (oo) in s s {H®°°), where all but q copies of H have been 
traced out. Infinite direct products of Hilbcrt spaces, as 
well as of states, are mathematically well-defined. The 
Hilbert space 7i®°° is the infinite direct product of H, 
defined as an inductive limit [17, 18]. The infinite direct 
product of a state p, p®°° , is the unique state on H.® 00 
such that if Aj = 1 for all but a finite number m of Ai, 



To show that A max (P g +„(j4 ® 1®")) is non-increasing 
with n, we use the convexity of A max over B H : 

A max (P g+ „+i(^®l® ( " +1) )) 
= A max (P g+n+1 (P g+n (A®ll®")®l) 



< 



, T ~~TTu E A max (pt(P 9+ „ (j 4®l®")®l)P OT ) 



= A max (P g+ „(A®l®")®l) 
= A max (P g+ „(A®l®")). 

This proves that the convergence of the right-hand side 
of eq. (2) is monotonically decreasing, so that the right- 
hand side of (2) is an upper bound on the left-hand side 
of (2) for all finite values of n. 

From the above calculations, we see that the optimal 
state p in max p Tr^p® 9 ] is approximated by 

p%\ « Tr„[P g+ „(|V>„)<V„|)] - Tr„ |V„)<Vn|, 
with V„ = Vw(P,+n(A ® 1®")). Therefore, 
p pt = lim Tr g+n _i \ip n )(ip n \. 



□ 



D. Simplifications for the Calculation of 



Tr[p®°° {® i A i )] = '[[Tr\pA j ]. 

As a result, we can express the states / p® q dp(p) as a 
partial trace of symmetric states in S S (H® 00 ): 

maxTr[Ap®«] 
p 



In this section we will show that in actual calculations 
of v q {<&) we do not have to deal with the (very large) 
space 

ft®(<3+»). Henceforth we will assume that A corre- 
sponds to v q {§), via the relation (1). Furthermore, for 
q > 2, we will assume that A is already the symmetric 
part of the originally obtained matrix. Thus, henceforth, 
A is symmetric and Hcrmitian. 
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The first thing we note is that the symmetrisation P q + n 
in P q + n (A ® 1®"), which in general involves a sum over 
all (q + n)\ permutations, can be replaced by a symmetri- 
sation over only C| + ™ permutations, provided A is sym- 
metric (or is made so, using P q ). Indeed, both factors A 
and 1®" already are symmetric, hence the only permu- 
tations necessary are those that redistribute the q indices 
of A among the q + n indices (considering the indices of 
A to be indistinguishable). 

The most important issue here is to simplify the cal- 
culation of the maximal eigenvalue of P q+n (A ® 1®"), 
and we can actually do that by exploiting the full per- 
mutation symmetry of that matrix. It is well-known that 
the eigenspaces of a fully symmetric matrix have to be 
themselves symmetry- invariant. Indeed, if A is symmet- 
ric then W e S n : P^AP^ = A. Inserting this in the 
eigenvalue equation Ax = Xx gives PlAP^x = Xx, or, 
using unitarity of P v , AP^x = XP^x. In other words, 
W e S n , P^x is in the same eigenspace as i. It fol- 
lows that every eigenspace is an S n -invariant subspace of 

The irreducible S^-invariant subspaccs correspond to 
the symmetry classes of S n . The best- known are the to- 
tally symmetric and the totally antisymmetric class, but, 
of course, there are many others, each one corresponding 
to a specific standard Young tableau of S n . The dimen- 
sion of an SVj-invariant subspace depends on the dimen- 
sion d of the underlying Hilbert space 7i. For the totally 
symmetric subspace, the dimension is C„ + " _1 , and for 
the totally antisymmetric subspace it is C^, which is zero 
if n > d. In fact, the dimension of any invariant subspace 
corresponding to a Young tableau of height larger than d 
is zero, so we only need to consider Young tableaux with 
height at most d. Still, the total number of invariant 
subspaces for a given n grows exponentially with n. 

If U is a unitary matrix (depending on d and n only) 
whose columns are basis vectors for the irreducible S n - 
invariant subspaces, then for a symmetric A, U'AU is 
block-diagonal, every block corresponding to one of the 
invariant subspaces. Likewise, for general A, P n {A) is 
unitarily equivalent with the matrix consisting of the di- 
rect sum of all the diagonal blocks of W AU . So, in order 
to calculate the maximal eigenvalue of a symmetric ma- 
trix, one needs only to calculate the maximal eigenvalues 
of each of the mentioned diagonal blocks, and then take 
the maximum of those. 

For every symmetry class, therefore, one can con- 
struct a matrix Uk (for the k-th symmetry class) whose 
columns are basis vectors of the corresponding invariant 
subspace, so that is the corresponding diagonal 

block. Moreover, if A is not symmetric, lllAUk is the cor- 
responding diagonal block for P n (A). Hence, in practice, 
the symmetrisation operation P„ need not be performed 
at all. This is good news, but the fact that the number 
of symmetry classes grows exponentially with n is still a 
nuisance. 

Quite fortunately, when we specialise to the case we 
are really interested in, namely the calculation of the 



maximal output purity of a channel, we only need to 
consider the totally symmetric subspace. Recall that 
i/f($) = max p Tr[$(p)«] = max. p Tr[Ap®«]. The func- 
tion Tr^p) 9 ] is convex, so its maximum over the state 
space will be achieved in a pure state p. Hence TrfAp® 9 ] 
will also be maximal for a pure p. On the other hand, 
Theorem 1 states that p opt = lim^oo Tr 9+ „_i \ip n )(i[> n \ 
with tp n = -0max(P 9 +™(^ <8> 1®")). The only way in 
which /9 opt can tend to a pure state, say V'opt, is if ip n 
tends to a tensor product ^ opt ® tp' (some -0'; this will of 
course have to be a tensor power of V> pt)- As a conse- 
quence, ip n must tend to a totally symmetric state. So, 
X ma ,x(S q +n(A ® l®™)S 9+ „) tends to the correct solution 
(where S 9 +„ is the projector on the totally symmetric 
subspace). Furthermore, it does so in a non-increasing 
way as well. Since S q+n is a projector, it is a contrac- 
tion. This, together with the fact that A is Hcrmitian 
here, is used for the inequality in: 

A max (S 9+ „ + i(A® 1®"® 1) 

~>q+n+l ) 

= X max (S q + n +i (S q + n (A ® l®™)S 9 +„ ® l)Sq+„+i) 

< X max (S q+n (A ® l®")S 9+ „ ® 1) 

= X max (S q + n (A Cg) l®")S 9 +„), 

and that is what we needed to show. We conclude that 
for the case of the MOP there is no point in considering 
other subspaces than the totally symmetric one. 

A final simplification is now possible. To reduce the 
dimension of the calculation we of course do not calcu- 
late S q+n (A ® l 0n )S 9+ „ directly, but rather P q+n {A ® 
l® n )P q+n , which is a matrix over a S(d,q + n)- 
dimensional Hilbert space H' (recall that P n is the ma- 
trix whose columns span the totally symmetric subspace 
S(7Y®"), while S q+n is the projector on it). There is 
no actual need to compute this over 

H®(q+n) directly, 

because this calculation is just a linear mapping from 
B(H® q ) to B(H'). 

We first introduce some additional notations. Recall 
that H is a d-dimensional Hilbert space and vectors in 
are indexed by (i) := . . . ,i n ). Every ij thus 
takes integer values from 1 to d. If x is a symmetric vector 
in then for any ir e S n , — x v (i)- This symme- 

try induces an equivalence relation on the set of indices, 
(i) ~ 7r(z) , and we will use as equivalence class represen- 
tative for an n-dimensional index (i) the <i-dimcnsional 
index [k] := [fci, . . . , kd] where kj is the number of times 
the value j occurs in the index (i). We denote this by 
[k] = #(«)■ Obviously, J2[k] := J2j =1 kj = n. One eas- 
ily sees that [k] can assume S(d, n) possible values. The 
size of the class corresponding to [k] is the value of the 
multinomial coefficient := n\/k\\ . . . fc^!. We adopt 
the common convention to take the value of the multino- 
mial coefficient to be whenever one or more of the kj 
is negative. 

Using these notations, an explicit form for P n is: 

{Pn){i),[k\ = S#(i),{k]C[ k ]- 
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The normalisation constant c^~] can be obtained from the 
requirement PjP n = 1: 

4] = #{(i):#(i) = [fc]} = q" fe] . 

Here #Z denotes the cardinality of the set Z. Thus 

{ P n)(i),[k] = S #(i),lk]( C [k])~ 1/2 - 

We will denote the quantity we want to calculate by 
Q n (A) :=pt +n (A®l®»)P, +n . 

Inserting the above expression for P q + n gives, with (i) 
and (j) now being (g + n)-dimensional indices, 

Qn(A) m i] = 6 #d),lk] 5 #(j),[l]( C [k\ nC [lt n y 1/2 

X ^(il l - l i(),l)'l l -j g ) ^»9+l,j 9 +l ■ • • $iq+n , jq+n 

(i').Ci') 

X E 5 #( i ' i ").[fe] <Wi")>[*]> 
(*") 

where (i') and (/) are (/-dimensional, (i") is n- 
dimensional, and (i'i") denotes the concatenation of (i') 
and (i"). Since #(ab) = #(a) + #(6), the sum over (i") 
reduces to 

S #(3'i")M = C [fc]-#(i') S lk]-#(i')M-#(j')- 

Thus, with [u] and [u] being d-dimensional indices whose 
sum equals q, 

(C$ n C$y/i Q n (A) mi] 

= S j4 (^),0'') c Tfe]-#( 4 ') *[fc]-#(i').w-#(j') 

(i').O'') 

KM 

x ]T {A (W) :#(i') = [«],#0") = M} 
(i').O') 

= S ^[fcl-M *[*]-[«]. M-M (^M^h) 172 
[ U ],M 

x (P}AP q ) [u]i[v] . 

This is final result. 

Note that the presence of the Kronecker delta in this 
formula implies that Q n (A) will have the same number of 
diagonals as P^AP qi which is 2S(d,q) — 1. The number 
of rows of Q n (A) is S(d, q + n), as we have noted before. 
If we treat Q n (A) as a sparse matrix, its storage require- 
ments are therefore just under (2S(d, q) — l)S(d, q + n). 
In function of n, this is of the order 0(n d_1 ). 



E. Numerical Study 

While our interest in the characterisation of v q {<&) by 
Theorem 1 is mainly theoretical — the hope is that it will 
ultimately lead to proving (or disproving) multiplicativ- 
ity of vi — it is also clear that Theorem 1 can be imme- 
diately converted to an algorithm for the global maximi- 
sation of Tr[Ap® 9 ] . Since this is a non-convex function of 
p it is most noteworthy that we have an algorithm here 
that is impervious to local maxima, unlike so many other 
optimisation algorithms. 

We have performed a preliminary numerical study on 
the convergence behaviour of the sequence fj, n (A), for the 
calculation of the maximal output 2-purity of qubit 
channels (d — 2 and q = 2). The matrix Q n (A) is in this 
case a banded matrix with 2 diagonals on both sides of 
the main diagonal. In all our experiments, the approx- 
imation error fJ, n (A) — fx^^A) was of the order 0(l/n). 
In fact, it turned out that the approximation error in 
function of n could itself be approximated to very good 
precision by a function of the form l/(an + 6), so much so 
that it proved possible to obtain a much better value of 
Moo (A) by extrapolating from a finite sequence of /i n (A) 
values. This extrapolation technique is similar in spirit to 
Aitkcn's method for accelerating the convergence of the 
power method for eigenvalue calculations. The obtained 
accuracy was frequently of the order of 10 -10 (depend- 
ing on the particular A), while the total running time 
was less than 1 second on a standard PC. More details 
on these experiments will be presented elsewhere. 

The time complexity of the algorithm is 0(n 2 ), as it 
is dominated by the calculation of the maximal eigen- 
value A max of a banded matrix of fixed width. We used 
the Matlab routine eigs (X, 1 , 'LR' ) (which calls the 
ARPACK library) to obtain the maximal eigenvalue of 
a Hermitian matrix X. 



III. CONCLUSION 

We have presented a global optimisation algorithm for 
a particular problem occurring in quantum information 
theory. The algorithm is based on representing the opti- 
misation problem as an eigenvalue problem over a larger 
space. The key to this representation is the Quantum 
de Finetti theorem, which is quite surprising given that 
its typical usage is in the study of the foundations of 
quantum mechanics. A lot of work remains to be done, 
both on the theoretical side (given our interest in exploit- 
ing Theorem 1 for settling the multiplicativity issues of 
v q ) and on the algorithmic one. It could very well be 
that Theorem 1 is just the simplest case of a larger class 
of optimisation problems that in a similar way can be 
transformed to an eigenvalue problem. We also believe 
this work will be of relevance to the study of higher-order 
tensors [5]. 
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